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Abstract 

Observations of light isotopes in cosmic rays provide 
valuable information on their origin and propagation in the 
Galaxy. Using the data collected by the AMS-01 experi- 
ment in the range '^0.2- l.SGeVnucleon"^, we compare 
the measurements on ^H, ^H, '^He, and ^He with calcu- 
lations for interstellar propagation and solar modulation. 
These data are described well by a diffusive-reacceleration 
model with parameters that match the B/C ratio data, indi- 
cating that He and heavier nuclei such as C-N-0 experience 
similar propagation histories. Close comparisons are made 
within the astrophysical constraints provided by the B/C ra- 
tio data and within the nuclear uncertainties arising from er- 
rors in the production cross section data. The astrophysi- 
cal uncertainties are expected to be dramatically reduced by 
the data upcoming from AMS-02, so that the nuclear uncer- 
tainties will likely represent the most serious limitation on 
the reliability of the model predictions. On the other hand, 
we find that secondary-to-secondary ratios such as ^H/'^He, 
®Li/^Li or i°B/^^B are barely sensitive to the key propaga- 
tion parameters and can represent a useful diagnostic test for 
the consistency of the calculations. 

Keywords cosmic rays — acceleration of particles — nu- 
clear reactions, nucleosynthesis, abundances 

1 Introduction 

Secondary Cosmic Ray (CR) isotopes such as ^H, '^He and 
Li-Be-B are believed to be produced as a results of of nu- 
clear interactions primary CRs such as ^H, ^He or C-N-0 
with the gas nuclei of the interstellar medium (ISM). The 
secondary CR abundances depend on the intensity of their 
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progenitors nuclei, their production rate and their transport 
in the turbulent magnetic fields (Strong et al. 2007). Sec- 
ondary to primary ratios such as ^H/'*He, •^He/'^He or B/C 
be used to study the CR propagation processes in the Galaxy. 
The B/C ratio is widely used to determine the key parame- 
ters of propagation models. In fact the B/C ratio is measured 
by several experiments between ~ lOOMeV and ~ 1 TeV of 
kinetic energy per nucleon. The CR propation physics is 
also connected with the indirect search of dark matter parti- 
cles. In this context the CR propagation models, once tuned 
to agree with the B/C ratio, are used to compute the sec- 
ondary production for other rare species such as p or d, 
that provides the astrophysical background for the search of 
new physics signals (Donatoetal. 2008; Evoli et al. 2011; 
Salati et al. 2010). Clearly, understanding the CR propa- 
gation processes is crucial for modeling both the CR sig- 
nal and the background. Furthermore, these studies as- 
sume that all the CR species experience the same propaga- 
tion effects in their journey thourghout the ISM (Putze et al. 
2010; Trottaetal. 2011). It is therefore important to test 
the CR propagation with nuclei of different mass-to-charge 
ratios. This issue of the universality of CR propagation his- 
tories was also studied in Webber (1997) and, recently, in 
Coste etal. (2011). 

In this work we use the recent AMS-01 observations for 
the ^H/^He and '^He/''He ratios and compare them with 
the expected ratios based on interstellar and heliospheric 
propagation calculations. The aim of this work is to deter- 
mine wether the AMS-01 observations are consistent with 
the propagation calculations derived from heavier nuclei 
(mainly from B/C data). This consistency is inspected 
within two classes of model uncertainties: the astrophysical 
uncertainties, which are related to the knowledge of the CR 
transport parameters given by the B/C ratio, and the nuclear 
uncertainties, which arise from the and ■^He production 
cross sections. 
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Fig. 1 Left: energy spectra of CR proton, deuteron (divided by 20), ^He and ''He. Right: isotopic ratios ^H/^'hie and ''He/'^He. 
Calculations are compared with the data from AlVIS-Ol (Aguilar et al. 2002, 2011), IMAX (Reimeretal. 1998; De Nolfo et al. 2000), 
SMILI (Ahlen et al. 2000), BESS (Wang et al. 2002), Hatano et al. (1995), Leech & O'Gallagher (1978), Webber & Yushak (1983). 



2 Observations 

The AMS-01 experiment operated successfully in the STS- 
91 mission on board the space shuttle Discovery. The 
spectrometer was composed of a cylindrical permanent 
magnet, a silicon micro-strip tracker, time-of-flight scin- 
tillator planes, an aerogel Cerenkov counter and anti- 
coincidence counters. The performance of AMS-01 is 
described elsewhere (Aguilar et al. 2002). Data collec- 
tion started on 1998 June 3 and lasted 10 days. AMS- 
01 observed cosmic rays at an altitude of ~380 km dur- 
ing a period, 1998 June, of relatively quiet solar activ- 
ity. Results on isotopic spectra have been recently pub- 
lished in (Aguilar et al. 2011) with the ratios ^H/*He, 
^HeZ-'He, *^Li/^Li, ^Be/CBe+^"Be) and ioB/"B in the 
range ^ 0.2 — 1.5 GeV of kinetic energy per nucleon. Fig. 1 
shows the AMS-01 energy spectra of proton, deuteron, 
helium isotopes, and the ratios ^H/^He and '^He/''He. 
The other data come from balloon borne experiments 
IMAX (Reimeretal. 1998; De Nolfo etal. 2000), SMILI 
(Ahlen et al. 2000), BESS (Wang et al. 2002), Hatano et al. 
(1995), Leech & O'Gallagher (1978) Webber & Yushak 
(1983). 

The AMS-01 observations are made in a period, 1998 
June, of relatively quiet solar activity, and the particle 
recorded are free from any atmospheric induced back- 
ground. Furthermore, the AMS-01 material thickness 
between the top of the payload and the active detector 
amounted to -^Sgcm"^ which is considerably less than 
that of previous balloon borne experiments (~ 9-20 g cm^^ 



of top-of-instrument material plus ~5gcm^^ of residual 
atmosphere). Also important, for the aims of this work, 
is to realize that high precision data are currently flowing 
in from two active projects PAMELA and AMS-02, both 
operating in space. In particular, the data forthcoming by 
are expected to provide a dramatic improvement in our un- 
derstanding of the CR transport processes and interactions 
(Tomassetti & Donato 2012; Coste et al. 201 1; Oliva 2008). 

3 CR Transport and Interactions 

Galactic CR nuclei are believed to be accelerated by particle 
diffuse shock acceleration mechanisms occurring in galac- 
tic sites such as supernova remnants (SNRs). Their prop- 
agation in the ISM is dominated by particle transport in 
the turbulent magnetic field and interactions with the mat- 
ter, that is generally described by a diffusion-transport equa- 
tion including source distribution functions, magnetic dif- 
fusion, energy losses, hadronic interactions, decays, dif- 
fusive reacceleration and convective transport (the latter 
is not considered in this work). Models of CR propa- 
gation in the Galaxy employ fully analytical (Thoudam 
2008; Tomassetti 2012), semi-analytical (Jones et al. 2001; 
Putze et al. 2010), or fully numerical calculation frame- 
works (Di Bernardo et al. 2010; Strong et al. 2007). The 
present work relies on the dijfusive-reacceleration model 
implemented with GALPROP-v50 . Ip', which numeri- 
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Fig. 2 Left: Cross section parametrizations (Cucinottaet al. 1993) for the channels "^He— ^'^He, "^He— j^'^H, ''He^^H andp+p — >• tt+^H 
(multiplied by 20). The data are encoded as ME72: Meyer (1972), WE90: Webber (1990), NI72: Nicholls et al. (1972), GL93: 
Glagolevetal. (1993), AB94: AbduUin et al. (1994), BLOl: Blinov et al. (2001), LM49: Lebowitz & Miller (1969), Right: B/C ra- 
tio from the CR propagation model of Table 1 and l-a uncertainty band. Data are from TRACER (Obermeier et al. 2011), AMS-01 
(Aguilar et al. 2010), CRIS (De Nolfo et al. 2003), CREAM (Ahn et al. 2008), HEAO (Engelmann et al. 1990), Orth et al. (1978), and 
Lezniak & Webber (1978). 



cally solves the cosmic ray propagation equation for a cylin- 
drical diffusive region with a realistic interestellar gas distri- 
bution and source distribution. 

3 . 1 Diffusive-Reacceleration Model 

The propagation equation of the diffusive-reacceleration 
model for a CR species j is given by: 



dt 



d_ 

dp 
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where Mj ~ dNj/dVdp is the CR density of the species i 
per unit of total momentum p. 

The source term, q^°^ = q^"^ + qj'^, includes the pri- 
mary acceleration spectrum {e.g. from SNRs) and the term 
arising from the secondary production in the ISM or decays. 
The primary spectrum is gj" = qo [R/ Rb)~'^ , is normal- 
ized to the abundances, 5°, at the reference rigidity Kb- The 
injection spectral indices are the same for all the primary el- 
ements. The source spatial distribution in the galactic disc 
is extracted from SNR observations. The secondary produc- 
tion term, q^'^ = Tlik-^f^^k^j' describe the products of 
decay and spallation of heavier CR progenitors with numer 
density 7\4- For collisions with the interstellar gas: 



Tk^j = I3kc 



at^AE,E')dE' , 



(2) 



where riism are the number densities of the ISM nuclei, 
nn ~ 0.9 cm~'^ and uhc ~ 0.1 cm~^, and cF^^™j are the 
fragmentation cross sections for the production of a j-type 
species at energy E from a fc-type progenitor of energy E' in 
H or He targets. = jijC {nH(jf}j + nHeCT]%) + ^ 
is the total destruction rate for inelastic collisions (cross sec- 
tion (7*°*) and/or decay for unstable particles (lifetime r). 
The spatial diffusion coefficient D ~ D{r,p) is taken as 
spatially homogeneus and rigidity dependent as D{R) = 
/3Do (R/Rq)^ , where R = pc/Ze is the rigidity, Dq fixes 
the normalization at the reference rigidity Rq, and the pa- 
rameter S specifies its rigidity dependence. Diffusive reac- 
celeration is described as diffusion process acting in mo- 
mentum space. It is determined by the coefficient Dpp for 
the momentum space diffusion: 



^p^VA 
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where va is the Alfven speed of plasma waves moving 
through the ISM. The last term describes Coulomb and ion- 
ization losses by means of the momentum loss rate pj. 

GALPROP solves the steady-stady equation dAf^/dt = 
for a cylindrical diffusion region of radius r,nax and half- 
thickness L with boundary conditions A/}(r=rma^)=0 and 
N'j{z=±L)=0. The local interstellar spectrum (LIS) is then 
computed for each species at the solar system coordinates 

= 8.5 kpc and z = 0: 



LIS 



(i?)==^AG-(r0,p), 



(4) 
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Fig. 3 Astrophysical (left) and nuclear (right) uncertainty bands for the predicted ratios ^H/''He, '^He/'^He and ^H/^He in compar- 
ison with the AlVIS-Ol data. Other data are from IMAX (Reimer et al. 1998; De Nolfo et al. 2000), SMILI (Ahlen et al. 2000), BESS 
(Wang et al. 2002) Hatano et al. (1995), and Webber & Yushak (1983). 



Table 1 Propagation parameter set. 



Parameter 


Name 


Value 


Injection, break value 


Rb [GV] 


9 


Injection, index below Rb 




1.82 


Injection, index above Rb 




2.36 


Diffusion, magnitude 


Do [cm^ s"^] 


5.75 • 10 


Diffusion, index 


5 


0.34 


Diffusion, ref. rigidity 


Ro [GV] 


4 


Reacceleration, Alfven speed 


VA [km s^"'^] 


36 


Galactic halo, radius 


7^ [kpc] 


20 


Galactic halo, height 


Zh [kpc] 


4 


Solar modulation parameter 


(p [MV] 


500 



where A is the mass number and the flux is given in 
units of kinetic energy per nucleon E. For the descritpion we 
adopt the "conventional model" which finely reproduces the 
CR elemental fluxes at intermediate energies of ~ lOOMeV 
- 100 GeV per nucleon. 

3.2 Heliospheric Propagation 

CRs in the solar neighbourhood undergo convection, dif- 
fusion and energy changes as results of the expansion of 
the solar wind. To describe the solar modulation effect, 
we adopt the so-called force-field approximation that arise 
from the a spherically symmetric description of the helio- 
sphere (Gleeson & Axford 1968). The correspondence be- 
tween the (modulated) top-of-atmosphere spectrum, <I)'^0^, 
and the (unmodulated) LIS spectrum of 4, is expressed 
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by the effective parameter </> (GV) through; 




where the LIS and TOA energies per nucleon are related by 
E^^ = E^'~^^ + The main parameters of the model are 
listed in Table 1. The remaining specifications are as in the 
file galdef_50p_599278 provided with the package. 

3.3 Fragmentation Cross Sections 

To compute the secondary nuclei production rate from the 
disintegration of the heavier CR nuclei, a large amount of 
cross section estimates is required. The accuracy of the cal- 
culated secondary spectra therefore depends on the relia- 
bility of the production and destruction cross sections em- 
ployed. 

The production of and '^He isotopes is mainly due 
to collision of ^He nuclei. The "^He isotopes are also pro- 
duced via decay of trithium ('^H ^> '^He) which, in turn, is 
predominantly created by "'He spallation. The most relevant 
projectile -^fragment processes for the ''He abundance are 
^He — > ■'He and "'He — ■'H. The main deuteron production 
channel is "'He — > ^H. Low energy deuterons are also created 
by the fusion reaction p + p tt+^H acting between ^ 300 
and ~900MeV of the proton energy (Meyer 1972). Al- 
though the p-p fusion cross section is very small (cr 3 mb), 
this reaction contributes appreciably to the -^H abundance 
because of the large CR proton flux. Spallation of heavier 
nuclei (C, O, Fe) give a minor contribution, roughly < 10 % 
of their fractional abundance. For all these channels the to- 
tal inclusive reaction can be realized in a number of possible 
final states. We employ the phenomenological parametriza- 
tion of Cucinotta et al. (1993). These parametrizations are 
shown in Fig. 2 together with the accelerator data. For ^H 
and ■'He, the partial contributions of break-up (B) and strip- 
ping (S) reactions are shown separately. 

For all these processes we have assumed that the frag- 
ment is ejected with the same the kinetic energy per nucleon 
E of the projectile, E'. This straight-ahead approximation, 
expressed by a{E, E') « a{E)S{E — £"), has been vaH- 
dated within some percent of accuracy for reactions involv- 
ing Z > 2 nuclei (Kneller et al. 2003) and for lighter species 
(Z < 3) (Cucinotta et al. 1993). For the p-p fusion channel, 
the kinetic energy per nucleon is not conserved due to the 
kinematic of the reaction: the energy of ejected deuterons 
is sistematically lower than that of the proton of a factor 
~ 4 on average. Using a straight-ahead fashion, we write 
a{E,E') « a{E)5{E - £,E'), where ^ =4 is flie aver- 
age inelasticity of the ^H. The p-p fusion contributes to the 
^H flux at energies below 250 MeV nucleon^ For nu- 
clear reactions involving heavier (Z > 2) nuclei, we use the 



default cross section parametrization of GALPROP. To ac- 
count for CR collisions with the interstellar helium 10% 
of the ISM) the parametrization of Ferrando et al. (1988) 
is used. Calculation of CR spectra and ratios ^H/"'He and 
^He/'^He are shown in Fig. 1 for the modulation intensity of 
(/) = 500MV. 

4 Model Uncertainties 

We consider two classes of uncertainty in the model esti- 
mates. The astrophysical uncertainties are those associated 
with the transport parameters constrained by the B/C ra- 
tio. The relevant parameters for the secondary productions 
are S, va and the ratio Dq/L. We perform a grid scan in 
the parameter space {S,va,Dq/L} by running GALPROP 
multiple times. The other parameters, e.g. the source pa- 
rameters and the modulation potential, are kept fixed. For 
each parameter configuration, we select the B/C-compatible 
models within one sigma of uncertainty in the statis- 
tics. We use B/C data from HEAO (Engelmann et al. 1990), 
CREAM (Ahnetal. 2008), AMS-01 (Aguilar et al. 2010) 
and Orth et al. (1978). Fig. 2 illustrates the uncertainty band 
derived by this procedure. Note that this method has severe 
limitations and allows to simultaneusly explore only some 
parameters. More robust strategies require advanced statis- 
tical tools, see e.g. Trotta et al. (2011), Putze et al. (2010) 
and in particular Coste et al. (201 1). However the purpose 
in this work is estimating the parameter uncertainties rather 
than determining their exact values. 

The nuclear uncertainties on the ^H and ''He calculations 
are those arising from uncertainties in their production cross 
sections. In order to estimate the cross section uncertainties 
using the information from the measurements, we re-fit the 
parametrizations with the data to determine their overall nor- 
malizations and associated errors. The uncertainty bands are 
shown in Fig. 2 for the main reactions of ^H and ''He pro- 
duction. These uncertainties can be directly translated into 
error bands for the predicted ratios ^H/^He and ■'He/^He. 
These uncertainty bands are shown in Fig. 3. 

The AMS-01 data agree well with calculations within the 
astrophysical band, indicating consistency with the propa- 
gation picture arising from the B/C analysis. It is also clear 
that Z < 2 nuclear ratios carry valuable information on the 
transport parameters, i.e., they can be in principle used to 
tighten the constraints given by the B/C ratio. On the other 
hand, the nuclear uncertainties represent an intrinsic limi- 
tation on the accuracy of the model predictions, as reported 
in the right panels of the figure. Only precise cross sec- 
tion data or more refined calculations may pin down these 
uncertainties. Unaccounted errors or systematic biases in 
cross section estimates cause errors on the predicted ratios 
which, in turn, may lead to a mis-determination of the CR 
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transport parameters. Given the level of precision expected 
from PAMELA or AMS-02, systematic errors in the cross 
section data may represent the dominant source of uncer- 
tainty for the model predictions of light CR isotopes. A 
strategy to check the model consistency with CR data is the 
use of secondary to secondary ratios such as ^H/^He. In 
fact, the and ^He isotopes have the same astrophysical 
origin and simlar progenitors (mainly ^He). Thus, their ra- 
tio is almost unsensitive to the propagation physics and can 
be used to probe the net effect of the nuclear interactions. 
In fact, a mis-consistency between calculations and ^H/'^He 
data would indicate the presence of systematic biases in the 
cross sections that cannot be re-absorbed by the propagation 
parameters. As illustrated in the bottom panels of Fig. 3, the 
tight astrophysical uncertainty band (left) indicates a little 
discrepancy between data and model which can be under- 
stood if one consider the nuclear uncertainty (right), which 
is dominant for the ^H/'^He ratio. Similarly, the use of other 
ratios such as ^Li/^Li, Li/Be or ^"B/^^B can represent a use- 
ful diagnostic tool for testing the overall consistency of the 
model. 



5 Conclusions 

We have compared new observations of the ^H/'^He and 
^He/'^He ratios in CRs made by the AMS-01 experiments 
with standard calculations of secondary production in the 
ISM. These ratios are well described by propagation models 
consistent with B/C ratio under a diffusive-reacceleration 
scenario, suggesting the He and heavier nuclei such as C- 
N-0 experience similar propagation histories. The accuracy 
of the secondary CR calculations relies on the accuracy of 
the cross sections employed. Given the level of precision ex- 
pected from AMS-02, the errors in the cross section data will 
likely represent the dominant source of uncertainty for the 
model predictions of rare CR isotopes such as or •^He. 
Similar issues may concern for ^ ""Li, ^'^'^°Be, or ^°'^^B. 
These errors may be reduced using more refined calculations 
or more precise accelerator data. For example, future obser- 
vations may require the departure from the straight-ahead 
approximation which is generally assumed in the CR prop- 
agation studies for light nuclei. Possible consistency checks 
for propagation models can be done using the secondary to 
secondary ratios, which are less sensitive to the astrophys- 
ical aspects of the interstellar CR propagation. The use of 
ratios such as ^H/^He, ^LW^U or ^^B/^^B can represent a 
useful diagnostic test for the reliability of the calculations: 
any CR propagation model, once tuned on secondary to pri- 
mary ratios, must correctly reproduce the secondary to sec- 
ondary ratios as well. Another model limitation concerns 
the solar modulation effect. Any more refined modeling re- 
quires to leave the force field approximation, that may be 



too simple to finely reproduce the future data of different 
A/Z isotopes. Our understanding of the CR heliospheric 
propagation may be dramatically improved by the AMS-02 
long-term observations of different CR species. 
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